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We combine complementary datasets to constrain dark energy. Using standard Big Bang Nucleosynthesis and the 
observed abundances of primordial nuclides to put constraints on Qq at temperatures near T ~ lMeV, we find the 
strong constraint f2<g(MeV) < 0.045 at 2a c.L. Under the assumption of flatness, using results from Cosmic Microwave 
Background (CMB) anisotropy measurements, high redshift supernovae (SN-Ia) observations and data from local cluster 
abundances we put a new constraint on the equation of state parameter wq < —0.85 at 68% c.L. 



Introduction. The discovery that the universe's 
evolution may be dominated by an effective cosmo- 
logical constant Q, is one of the most remarkable 
cosmological findings of recent years. An exceptional 
opportunity is now opening up to decipher the na- 
ture of dark matter [|| , to test the veracity of theories 
and reconstruct the dark matters profile using a wide 
variety of observations over a broad redshift range. 

One candidate that could possibly explain the ob- 
servations is a dynamical scalar "quintessence" field. 
One of the strong aspects of quintessence theories 
is that they go some way explaining the fine-tuning 
problem, why the energy density producing the ac- 
celeration is ~ 10 _120 Mp ; . A vast range of "tracker" 
(see for example ||[llj) and "scaling" (for example 
01" 01) quintessence models exist which approach 
attractor solutions, giving the required energy den- 
sity, independent of initial conditions. The common 
characteristic of quintessence models is that their 
equations of state, wq = p/p, vary with time whilst a 
cosmological constant remains fixed at wq=a = — 1. 
Observationally distinguishing a time variation in 
the equation of state or finding wq different from 
— 1 will therefore be a success for the quintessential 
scenario. 

We will here discuss observational constraints on 
general quintessence models. 

In Q we used standard big bang nucleosynthe- 
sis and the observed abundances of primordial nu- 
clides to put constraints on the amplitude of the en- 
ergy density, SIq, at temperatures near T ~ lMeV. 
The inclusion of a scaling field increases the expan- 
sion rate of the universe, and changes the ratio of 
neutrons to protons at freeze-out and hence the pre- 
dicted abundances of light elements. 

In |Q we have then combined the latest observa- 
tions of the Cosmic Microwave Background (CMB) 
anisotropics provided by the Boomerang |lq ], DASI 
0] and Maxima 0] experiments and the informa- 
tion from Large Scale Structure (LSS) with the lu- 
minosity distance of high redshift supernovae (SN- 



Ia) to put constraints on the dark energy equation 
of state parameterized by a redshift independent 
quintessence-field pressure-to-density ratio wq. We 
also made use of the Hubble Space Telescope (HST) 
constraint on the Hubble parameter h = 0.72 ± 0.08 




We will briefly review the results obtained in those 
paper in the next sections. 

Early-Universe Constraints from BBN. In the last 
few years important experimental progress has been 
made in the measurement of light element primor- 
dial abundances. For the 4 He mass fraction, Ufj e , 
two marginally compatible measurements have been 
obtained from regression against zero metallicity in 
blue compact galaxies. A low value Ypj„ = 0.234 ± 
0.003 § and a high one U He = 0.244 ± 0.002 § 
give realistic bounds. We use the high value in our 
analysis; if one instead considered the low value, the 
bounds obtained would be even stronger. 

Observations in different quasar absorption line 
systems give a relative abundance of deuterium, 
which is critical in fixing the baryon fraction, of 
D/H = (3.0 ±0.4) • 10~ 5 @. 

In the standard BBN scenario, the primordial 
abundances are a function of the baryon density 
r\ ~ fit, h 2 only. In order to put constraints on the 
energy density of a primordial field aT~ MeV, we 
modified the standard BBN code Q by including the 
quintessence energy component Qq. We then per- 
formed a likelihood analysis in the parameter space 
(Qbh 2 ,fiQ BN ) using the observed abundances Yg e 
and D/H. In Fig. [j] we plot the 1,2 and 3<r likeli- 
hood contours in the (Q^h 2 ,£Iq BN ) plane. 

Our main result is that the experimental data for 
4 He and D does not favour the presence of a dark 
energy component, providing the strong constraint 
^Q(MeV) < 0.045 at 2a (corresponding to A > 9 for 
the exponential potential scenario), strengthening 
significantly the previous limit of |8|;[l5| Qq (MeV) < 
0.2. The reason for the difference is due to the 
improvement in the measurements of the observed 
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Figure 1. 1,2 and 3er likelihood contours in the 
(f2f,/i 2 , fi<g(lMeV)) parameter space derived from 
4 He and £) abundances. 



abundances, especially for the deuterium, which now 
corresponds to approximately AiVgff < 0.2 — 0.3 ad- 
ditional effective neutrinos (see, e.g. whereas 
Rcf. p|,|l5|| used the conservative value AiV e g- < 1.5. 
One could worry about the effect of any underesti- 
mated systematic errors, and we therefore multiplied 
the error-bars of the observed abundances by a fac- 
tor of 2. Even taking this into account, there is still 
a strong constraint fiQ(MeV) < 0.09 (A > 6.5) at 
2a. 

Constraints on the Dark energy equation of state. 

The importance of combining different data sets in 
order to obtain reliable constraints on wq has been 
stressed by many authors (see e.g. Q, |2"H| , [p2| ) , 
since each dataset suffers from degeneracies between 
the various cosmological parameters and wq . Even 
if one restricts consideration to flat universes and 
to a value of wq constant in time then the SN-la 
luminosity distance and position of the first CMB 
peak are highly degenerate in wq and Oq , the energy 
density in quintessence. 

The effects of varying wq on the angular power 
spectrum of the CMB anisotropics can be reduced to 
just two . Firstly, since the inclusion of quintessence 
changes the overall content of matter and energy, 
the angular diameter distance of the acoustic horizon 
size at recombination will be altered. In flat models 
(i.e. where the energy density in matter is equal 
to VI m = 1 — ^q), this creates a shift in the peaks 
positions of the angular spectrum as 
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y = 
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+n Q (i + z) 3( - 1+w ^]- 1/2 dz. 

It is important to note that the effect is completely 



degenerate in the interplay between wq and CIq. 
Furthermore, it does not add any new features be- 
yond those produced by the presence of a cosmolog- 
ical constant |25[] , and it is not particularly sensitive 
to further time dependencies of wq. 
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Figure 2. CMB power spectra and the angular diam- 
eter distance degeneracy. The models are computed 
assuming flatness, Qfc = 1 — £Im — = 0). The 
Integrated Sachs Wolfe effect on large angular scale 
slightly breaks the degeneracy. The degeneracy can 
be broken with a strong prior on h, in this paper we 
use the results from the HST. 



Secondly, the time-varying Newtonian potential 
after decoupling will produce anisotropies at large 
angular scales through the Integrated Sachs- Wolfe 
(ISW) effect. The curve in the CMB angular spec- 
trum on large angular scales depends not only on the 
value of wq but also its variation with redshift. How- 
ever, this effect will be difficult to disentangle from 
the same effect generated by a cosmological constant, 
especially in view of the affect of cosmic variance 
and/or gravity waves on the large scale anisotropies. 

In order to emphasize the importance of degen- 
eracies between all these parameters while analyzing 
the CMB data, we plot in Figure 2 some degener- 
ate spectra, obtained keeping the physical density in 
matter J7m^ 2 , the physical density in baryons Vl^ti 2 
and 1Z fixed. As we can see from the plot, models de- 
generate in wq can easily be constructed. However 
the combination of CMB data with other different 
datasets can break the mentioned degeneracies. 

Type la supernovae, in particular, can be ex- 
tremely useful in this. Evidence that the universe's 
expansion rate was accelerating was first provided by 
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Figure 3. Contours of constant 1Z (CMB) and 
SN — la luminosity distance in the wq-Hm plane. 
The degeneracy between the two distance measures 
can be broken by combining the two sets of com- 
plementary information. The luminosity distance is 
chosen to be equal to di at z = 1 for a fiducial model 
with Q\ = 0.7, fi M = 0.3,ft = 0.65. (We note that as 
£Iq = 1 — VLm goes to zero the dependence of 1Z and 
d,L upon wq also become zero, as there is no dark 
energy present.) Figure taken from p4| 



two groups, the SCP and High-Z Search Team(|l|) 
using type la supernovae (SN-Ia) to probe the nearby 
expansion dynamics. SN-Ia seem to be good stan- 
dard candles, as they exhibit a strong phenomcno- 
logical correlation between the decline rate and peak 
magnitude of the luminosity. The observed appar- 
ent bolometric luminosity is related to the lumi- 
nosity distance, measured in Mpc, by = M + 
hlogdL(z) + 25, where M is the absolute bolomet- 
ric magnitude. The luminosity distance is sensi- 
tive to the cosmological evolution through an inte- 
gral dependence on the Hubble factor di = (1 + 
z) L (dz'/H(z',QQ,Wq), and can therefore be used 
to constrain the scalar equation of state. As can be 
seen in Figure ^ there is an inherent degeneracy in 
the luminosity distance in the Qm / wq plane; one can 
see that little can be said about the equation of state 
from luminosity distance data alone. However, the 
degeneracies of CMB and SNla data complement one 
another so that together they offer a more powerful 
approach for constraining wq. 

Table I shows the 1-a constraints on wq for dif- 
ferent combinations of priors, obtained in E4] af- 



Table 1 

Constraints on wq and Qm = 1 — using differ- 
ent priors and datasets. We always assume flatness 
and to > 10 Gyr. The la limits are found from 
the 16% and 84% integrals of the marginalized like- 
lihood. The HST prior is h = 0.72 ± 0.08 while 
for the BBN prior we use the conservative bound 
n b h 2 = 0.020 ±0.005. 



CMB+HST wq < -0.62 

0.15 < n M < 0.45 
CMB+HST+BBN -0.95 <iu Q < -0.62 

0.15 < Q M < 0.42 
CMB+HST+SN-Ia -0.94 < w Q < -0.74 

0.16 < n M < 0.34 
CMB+HST+SN-Ia+LSS w Q < -0.85 

0.28 < Q M < 0.43 



ter marginalizing over all remaining nuisance pa- 
rameters. The analysis is restricted to flat uni- 
verses. One can see that wq is poorly constrained 
from CMB data alone, even when the strong HST 
prior on the Hubble parameter, h = 0.72 ± 0.08, 
is assumed. Adding a Big Bang Nucleosynthesis 
prior, Q b h 2 = 0.020 ± 0.005, has small effect on 
the CMB+HST result. Adding SN-Ia breaks the 
CMB Qq — wq degeneracy and improves the upper 
limit on wq, with wq < —0.74. Finally, including 
information from local cluster abundances through 
(78 = (0.55 ± 0.1)r^ f ' 5 , where as is the rms mass 
fluctuation in spheres of Sh" 1 Mpc, further breaks 
the quintessential-degeneracy, giving wq < —0.85 at 
1-a. Also reported in Table I, are the constraints on 
Vtjw As we can see, the combined data suggests the 
presence of dark energy with high significance, even 
in the case CMB+HST. It is interesting to project 
our likelihood in the Qq — wq plane. Proceeding as in 
[ p7| , we attribute a likelihood to a point in the (CIm, 
wq ) plane by finding the remaining parameters that 
maximize it. We then define our 68%, 95% and 99% 
contours to be where the likelihood falls to 0.32, 0.05 
and 0.01 of its peak value, as would be the case for a 
two dimensional multivariate Gaussian. In Figure 4 
we plot likelihood contours in the (£Im, w q) plane for 
the joint analyses of CMB+SN-Ia+HST+LSS data 
together with the contours from the SN-Ia dataset 
only. As we can see, the combination of the data 
breaks the luminosity distance degeneracy. 

Conclusions. 
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Figure 4. The likelihood contours in the (£Im, wq) 
plane, with the remaining parameters taking their 
best-fitting values for the joint CMB+SN-Ia+LSS 
analysis described in the text. The contours corre- 
spond to 0.32, 0.05 and 0.01 of the peak value of the 
likelihood, which are the 68%, 95% and 99% confi- 
dence levels respectively. 



We have provided new constraints on the dark en- 
ergy by combining different cosmological data. We 
have examined BBN abundances in a cosmological 
scenario with a scaling field. We have quantitatively 
discussed how large values of the fractional density 
in the scaling field VIq at T ~ lMeV can be in agree- 
ment with the observed values of 4 He and D, assum- 
ing standard Big Bang Nucleosynthesis. The 2a limit 
f2<3(lMeV) < 0.045 severely constrains a wide class 
of quintessential scenarios, like those based on an ex- 
ponential potential. For example, for the pure expo- 
nential potential the total energy today is restricted 
to n Q = |n Q (lMeV) < 0.04. This result put strong 
constraints on the presence of quintessence during 
recombination. 

We have then provided constraints on the equa- 
tion of state parameter wq. The new CMB results 
provided by Boomerang and DASI improve the con- 
straints from previous and similar analysis (see e.g., 
||,||), with w Q < -0.85 at 68% c.l. (w Q < -0.76 
at 95% c.l.). We have also demonstrated how the 
combination of CMB data with other datasets is 
crucial in order to break the VIq — wq degeneracy. 
The constraints from each single datasets are, as ex- 
pected, quite broad but compatible with each other, 
providing an important consistency test. When com- 
parison is possible (i.e. restricting to similar priors 



and datasets), our analysis is compatible with other 
recent analysis on wq ( |30tl ). Our final result is per- 
fectly in agreement with the wq = — 1 cosmological 
constant case and gives no support to a quintessen- 
tial field scenario with wq > — 1. A frustrated net- 
work of domain walls or a purely exponential scaling 
field are excluded at high significance. In addition 
a number of quintessential models are highly disfa- 
vored, e.g. power law potentials with p > 1 and the 
oscillatory potential, to name a few. 

It will be the duty of higher rcdshift datasets, for 
example from clustering observations [p9| to point to 
a variation in w that might place quintessence in a 
more favorable light. 

The result obtained here, however, could be 
plagued by some of the theoretical assumptions 
we made. The CMB and LSS constraints can be 
weakened by the inclusion of additional relativis- 
tic degrees of freedom Q, by a background of 
gravity waves, of isocurvature perturbations or by 
adding features in the primordial perturbation spec- 
tra. These modifications are not expected in the 
most basic and simplified inflationary scenario but 
they are still compatible with the present data. The 
SN-Ia result has been obtained under the assump- 
tion of a constant- with-time wq > — 1. Inclusion of 
the region with wq < — 1 could affect our constraints 
( |32| ) . In |^4| we have shown that in general w e f t is a 
rather good approximation for dynamical quintessen- 
tial models since the luminosity distance depends on 
wq through a multiple integral that smears its red- 
shift dependence, and our results are therefore valid 
for a wide class of quintessential models. This 'numb- 
ing' of sensitivity to wq , first noticed by |3l[] , implies 
that maybe an effective equation of state is the most 
tangible parameter able to be extracted from super- 
novae. However with the promise of large data sets 
from Planck and SNAP satellites, opportunities may 
yet still be open to reconstruct a time varying equa- 
tion of state ]22] ]. 
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